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1 Introduction 

This paper is an expanded version of |2H] and [301 where there is presented an 
introduction to a point of view for discrete foundations of physics. In taking 
a discrete stance, we find that the initial expression of physical observation 
naturally occurs in a context of non-commutative algebra. In this way a 
formalism similar to quantum mechanics occurs first, but not necessarily with 
the usual interpretations. By following this line we show how the outlines of 
the well-known forms of physical theory arise first in non-commutative form. 
The exact relation of commutative and non-commutative theories raises a 
host of problems. 

The starting point for this investigation is the representation of calculus 
in a non-commutative framework. In such a framework derivatives are repre- 
sented by commutators, or more generally by products that satisfy the Jacobi 
identity and the Leibniz rule. If we take commutators [A, B] = AB — BA 
in an abstract algebra and define DA = [A, J] for a fixed element J, then 
D acts like a derivative in the sense that D(AB) = D(A)B + AD(B) (the 
Leibniz rule). As soon as we have calculus in such a framework, concepts of 
geometry are immediately available. For example, if we have two derivatives 
VjA = [A, J] and VkA = [A,K], then we can consider the commutator of 



these derivatives [Vj, V K ]A = V jV K A - V K V jA = [[J,K],A]. The non- 
commutation of derivations corresponds to curvature in geometry, and indeed 
we shall see that the earliest emergence of curvature in this context is the 
formal analog of the curvature of a gauge connection! 

For multivariable calculus we need variables Xi,X2, ■ ■ ■ , X n and elements 
Pi, P2, P n such that diA = dA/dXi = [A, Pi]. For a simplest representa- 
tion we shall assume the the commute with one another, and that the Pj 
commute with one another. Since we want diXj = 5ij (the Kronecker delta 
5ij is equal to one if i and j are equal and is zero otherwise), we must have 
the commutator equation [Xj,Pj] = 5^. Thus multivariable calculus in this 
non-commutative representation demands the commutation relations 

[Xi,Xj] = 

[^,^■1=0 
[X u Pj] = Sij 

These equations are the "flat background" for our non-commutative calcu- 
lus. The reader will note that this flat background has the same pattern of 
commutation relations as a bare form of quantum mechanics when the X 
variables are interpreted as position and the P variables are interpreted as 
momenta. In a certain sense this means that our considerations start in the 
quantum domain. Note that flat is a correct adjective, since the derivatives 
di all commute with one another. 

Let Ai be a collection of elements of this algebra. Define "covariant 
derivatives" with A; = Pi — Ai by the formula 

ViZ= [Z,Ai] =diZ-[Z,Ai]. 

Computing the curvature, one finds 

[Vi,Vj]Z=[[Ki,Kj],Z] 

and 

[A h Aj] = diAj - djAi + [A h Aj]. 

The reader will recognize this last expression as the formula for the curvature 
of a gauge connection. 
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In interfacing this formalism with physics we adopt the coupling equation 

dXi/dt = Ai = Pi - Ai. 

The reader will recognise this as the minimal coupling postulate in the con- 
text of Poisson brackets. Here we take it in the context of commutators or 
Poisson brackets, or a more general product satisfying the Jacobi identity and 
the Leibniz rule as described above. One retrieves the physics of a gauge field 
in this formalism. This is the essence of the pattern behind the Feynman- 
Dyson derivation of electromagnetism from commutation relations |2*H |2*T] , 
and its import is more general. Because the brackets can be interpreted 
as commutators or as Poisson brackets with special structure, the formalism 
can be seen in a multiplicity of contexts. Deeper relationships with curvature 
and metric are related to this shifting of contexts clS 3X6 relationships with 
quantum mechanics where the quantum formalism is obtained by the Dirac 
prescription of replacing Poisson brackets by commutators. We will discuss 
these issues in Section 5 of this paper. The organization of the paper is as 
follows. 

Section 2 of this paper we discuss the properties of the non- commutative 
discrete calculus that underlies our work. Here we begin with the considera- 
tion of a temporal operator J with the property that YJ = JY' for a "time 
series" X, X\ X", • • • . Thus XJ = JX', X'J = JX", This formalism 
for time series gives rise to the time derivative DA = [A, J] = AJ — J A = 
J A — J A = J (A! — A), a commutator representing a discrete derivative. 
Note that DA satisfies the Leibniz rule, a privilege not shared by the usual 
commutative discrete derivative. This section discusses the discrete ordered 
calculus (DOC) that arises from this idea and applies these ideas to a num- 
ber of situations. In particular, we consider the one variable case of the 
commutator equation [X, DX] = Jk and show that it leads to a Brownian 
walk, and that if we take the size of the time step into account, then the 
diffusion constant for a Brownian process arises naturally as k/2. We com- 
pare this with the usual derivation of the diffusion constant and the diffusion 
differential equation. We then compare this situation with the one dimen- 
sional Schrodinger equation, modeling it in relation to a diffusion process 
with complex amplitudes. In this viewpoint one sees that the step length of 
the diffusion process is the Compton wavelength associated with the mass 
for the particle, and the time is the Compton time. For the Planck mass this 
gives a step equal to the Planck length and a time interval equal to the Planck 
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time. We speculate on the relationship of this result to joint work with Pierre 
Noyes and others jHZj. We consider other time series that can be regarded as 
solutions to this Heisenberg relation, the problem of using more variables and 
a model that is related to a discrete version of the Feynman-Dyson derivation 
of electromagnetic formalism. 

Section 3 examines the consequences for a particle whose position - mo- 
mentum commutator is equated to a metric field. Here we see how the 
Levi-Civita connection (and implicitly differential geometric structure) comes 
naturally from the non- commutative calculus. This is a very general result 
and in section 4 we discuss it in a more axiomatic context as described in this 
introduction. This section discusses the intimate relationship between that 
Levi-Civita connection and the Jacobi and Leibniz identities that is revealed 
by our non- commutative calculus. In section 5 our stance leads to an inver- 
sion of the usual Dirac maxim "replace Poisson brackets with commutators" . 
If we replace commutators with Poisson brackets that obey a Leibniz rule 
satisfied by the commutators, then the dynamical variables will obey Hamil- 
ton's equations. Thus we can take Hamilton's equations as a classicization of 
our theory. Among other things, this point of view explains the appearance 
of the Levi-Civita connection in the abstract formalism. Interpreting with 
Poisson brackets, we obtain a new proof (via Jacobi identity) of the classical 
result that a Newtonian particle moving in generalized coordinates according 
to Lagrange's equations describes a geodesic in the Levi-Civita connection. 
Section 6 discusses the relationship of the discrete ordered calculus with q- 
deformations and quantum groups. We show that in a quantum group with 
a special grouplike element representing the square of the antipode, there is a 
representation of the discrete ordered calculus. In this calculus on a quantum 
group the square of the antipode can represent one tick of the clock. Then 
follows section 7 on networks and discrete spacetime. This section is an expo- 
sition of ideas related to spin networks and topological quantum field theory. 
As an early example we discuss the discretization of the Dirac equation in 
1 + 1 dimensional spacetime. It is our speculation that the approaches to 
discrete physics inherent in discrete calculus and in topological field theory 
are deeply interrelated. At the end of this section we outline this relationship 
in the case of a model for quantum gravity due to Louis Crane. Section 8 
is an appendix on the iterant approach to matrix algebra. We include this 
appendix to show how one can conceptualize matrix algebra from point of 
view of the discrete. Section 9 is a philosophical appendix discussing the 
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nature of foundations in mathematics and in physics. 

Remark. The following references in relation to non-commutative calculus 
are useful in comparing with our approach (TUTUIIISIES]- Much of the present 
work is the fruit of a long series of discussions with Pierre Noyes, and we will 
be preparing collaborative papers on it. The present paper is a summary 
for the proceedings of the ANPA Conference held in Cambridge, England 
in the summer of 2002. I particularly thank Eddie Oshins for pointing out 
the relevance of minimal coupling. The paper [33] also works with minimal 
coupling for the Feynman-Dyson derivation. The first remark about the 
minimal coupling occurs in the original paper by Dyson [3], in the context 
of Poisson brackets. The paper [TJJ is worth reading as a companion to 
Dyson. In the present paper we generalize the minimal coupling to contexts 
including both commutators and Poisson brackets. The reader can see the 
full generality of our approach by first reading this introduction and then 
going directly to sections 4 and 5. It is the purpose of this paper to indicate 
how non-commutative calculus can be used in foundations. 

Acknowledgement. Most of this effort was sponsored by the Defense Ad- 
vanced Research Projects Agency (DARPA) and Air Force Research Labo- 
ratory, Air Force Materiel Command, USAF, under agreement F30602-01-2- 
05022. Some of this effort was also sponsored by the National Institute for 
Standards and Technology (NIST). The U.S. Government is authorized to 
reproduce and distribute reprints for Government purposes notwithstanding 
any copyright annotations thereon. The views and conclusions contained 
herein are those of the author and should not be interpreted as necessarily 
representing the official policies or endorsements, either expressed or implied, 
of the Defense Advanced Research Projects Agency, the Air Force Research 
Laboratory, or the U.S. Government. (Copyright 2003.) It gives the author 
pleasure to thank Pierre Noyes, Clive Kilmister, Ted Bastin, Tony Deakin, 
Eddie Oshins, Basil Hiley, Keith Bowden, Arleta Giffor, Ashok Gangadean, 
Lynnclaire Dennis, Louis Licht and Sam Lomonaco for many conversations 
during the course of this work, and the Theory Group of the Stanford Linear 
Accelerator Laboratory for hospitality during the preparation of parts of the 
present paper. 
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2 Discrete Ordered Calculus 



In this section we recall the construction of an ordered version of the calcu- 
lus of finite differences DOC [21], [2E|- I n this calculus the Leibniz rule is 
satisfied, and so the calculus can be used in a variety of applications. 

In the abstract framework of this calculus, there are variables X, each of 
which connotes a time series 

x,x',x",.... 

Discrete unit time steps are indicated by the primes appended to the X. 
A general point in the time series at time t will be denoted by X t . By 
convention let the time step between successive points in the series be equal 
to 1 : 

At = 1. 

Then we can define the velocity at time t by the formula: 

v (t) = X t+l - X 1 . 

More generally, if X denotes position at a given time, then X' — X denotes 
the velocity at that time, where the phrase "at that time" must involve the 
next time as well. In a discrete context there is no notion of instantaneous 
velocity. 

Measure position, and you find X. Then measure velocity, and you get 
X' — X. Now measure position, and you get X' because the time has shifted 
to the next time in order to allow the velocity measurement. In order to 
measure velocity the position is necessarily shifted to its value at the next 
time step. In this sense, position and velocity measurements cannot commute 
in a discrete framework. This is the key physical idea that motivates our 
constructions. It was this idea, told to the author by Pierre Noyes, that led 
to our papers and particularly to [21] ■ 

The simplest interpretation of the variable X is that the time series values 
are numerical values, commuting with one another and with any operators 
that might be present in the associated mathematics or physics. In fact, we 
will often deal with situations where the X and the elements of the time 
series are in fact operators, not necessarily commuting with one another. At 
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the very least we will construct an algebra that mirrors the discrete non- 
commutativity of the operations of position and velocity measurement. 

Our project is to take this basic noncommutativity at face value and 
follow out its consequences. To this end we will formulate a calculus of 
finite differences that takes the order of observations into account. This 
formalization is explained below. 

To see most clearly the non-commutativity that is at the base of our 
considerations, let J denote the operation of shifting time by one increment. 
Thus we can envisage an algebra of operations that consists in commands 
like JX (measure X, then tick the clock). Note that we will agree to take 
the sequence of operations from right to left. Let \JX\ denote the "spatial 
evaluation" of this sequence of operations, obtained in general by performing 
all the instructions and then evaluating the spatial position. Thus 

\JX\ = X 

while 

\XJ\ = X' 

since when the clock ticks, the position shifts to the position at the next 
time. We see therefore, that XJ ^ JX. This is the first instance of non- 
commutativity in the physics of discrete space and time. From the point of 
view of spatial evaluation it is most convenient to declare the equation 

XJ = JX' 

since these two expressions have identical spatial evaluations. 

We can then define the DOC derivative by the equation 

DX = [X, J]= XJ - JX = JX' - JX = J(X' -X) = JdX 

where dX denotes the classical discrete derivative with unit time step. The 
key point about the DOC derivative is that it is a commutator, and conse- 
quently satisfies the Leibniz rule 

D(XY) = D(X)Y + XD(Y). 

This makes it possible to do discrete calculus in a way that is formally similar 
to classical calculus. We will repeat this structure more slowly now, first 
recalling the properties of classical discrete derivatives. 
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We begin by recalling the usual derivative in the calculus of finite differ- 
ences, generalised to a (possibly) non-commutative context. 

Definition. Let 

dX = X' - X 

define the finite difference derivative of a variable X whose successive values 
in discrete time are 

X,X',X",.... 

This dX is a classical derivative in the calculus of finite differences. It is 
still defined even if the quantities elements of the time series are in a non- 
commutative algebra. We shall assume that the values of the time series 
are in a possibly non-commutative ring R with unit. (Thus the values could 
be real numbers, complex numbers, matrices, linear operators on a Hilbert 
space, or elements of an appropriate abstract algebra.) This means that 
for every element A of the ring R there is a well-defined successor element 
A', the next term in the time series. It is convenient to assume that the 
ring itself has this temporal structure. In practice, one is concerned with a 
particular time series and not the structure of the entire ring. Moreover, we 
shall assume that the next-time operator distributes over both addition and 
multiplication in the sense that 

(A + B)' = A' + B' 

and 

(AB)' = MB'. 

An element c of the ring R is said to be a constant if d = c. 
Lemma 1. 

d(XY) = X'd{Y) + d(X)Y. 

Proof. 

d(XY) = X'Y' - XY 
= XV - X'Y + X'Y - XY 
= X'(Y' -Y) + (X' - X)Y 
= X'diY) + d{X)Y. 
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This formula is different from the usual formula in Newtonian calculus by 
the time shift of X to X' in the first term. We now correct this discrepancy in 
the calculus of finite differences by taking a new derivative D as an instruction 
to shift the time to the left of the operator D. That is, we take XD(Y) quite 
literally as an instruction to first find dY and then find the value of X. 
In order to find dY the clock must advance one notch. Therefore X has 
advanced to X' and we have that the evaluation of XD(Y) is 

X\Y' -Y). 



In order to keep track of this non-commutative time-shifting, we will write 

DX = J(X' - X) 
where the element J is a special time-shift operator satisfying 

ZJ = JZ' 

for any Z in the ring R. The time-shifter, J, acts to automatically evaluate 
expressions in the resulting non-commutative calculus of finite differences. 
We call this calculus DOC (for discrete ordered calculus). Note that J for- 
malizes the operational ordering inherent in our initial discussion of velocity 
and position measurements. An operator containing J causes a time shift in 
the variables or operators to the left of J in the sequence order. 

Formally, we extend the ring of values R (see the definition of d above) 
by adding a new symbol J with the property that AJ = J A' for every A 
in R. It is assumed that the extended ring R is associative and satisfies the 
distributive law so that J(A + B) = JA + JB and J(AB) = (JA)B for all 
A and B in the ring. We also assume that J itself is a constant in the sense 
that J' = J. 

The key result in DOC is the following adjusted difference formula: 
Lemma 2. 

D(XY) = XD(Y) + D(Y)X. 
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Proof. 

D(XY) 
= J(X'Y' - XY) 
= J(X'Y' - X'Y + X'Y - XY) 
= J(X'(Y' -Y) + (X' - X)Y 
= JX'(Y' -Y) + J(X' - X)Y 
= XJ(Y' -Y) + J(X' - X)Y 
= XD(Y) + D(X)Y. 

The upshot is that DOC behaves formally like infinitesimal calculus and 
can be used as a calculus in this version of discrete physics. In [21] Pierre 
Noyes and the author use this foundation to build a derivation of a non- 
commutative version of electromagnetism. Another version of this deriva- 
tion can be found in |27j. In both cases the derivation is a translation to 
this context of the well-known Feynman-Dyson derivation of electromagnetic 
formalism from commutation relations of position and velocity. 

Note that the definition of the derivative in DOC is actually a commu- 
tator: 

DX = J(X' -X) = JX' - JX = XJ - JX = [X, J}. 

The operator J can be regarded as a discretised time-evolution operator in 
the Heisenberg formulation of quantum mechanics. In fact we can write 
formally that 

X' = J~ X XJ 

since JX 1 = XJ (assuming for this interpretation that the operator J is 
invertible). Putting the time variable back into the equation, we get the 
evolution 

X t + At = j~l X tj 

This aspect can be compared to the formalism of Alain Connes' theory of 
non- commutative geometry jjj. 

In DOC, X and DX have no reason to commute: 
[X, DX] = XJ(X' - X) - J(X' - X)X = J{X'(X' -X)- (X 1 - X)X) 
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Hence 

[X, DX\ = J(X'X' - 2X'X + XX). 

This is non-zero even in the case where X and X' commute with one another. 
Consequently, we can consider physical laws in the form 

[X h DX j ]=g ij 

where is a function that is suitable to the given application. In [21] 
we show how the formalism of electromagnetism arises when g % i is 5^ , the 
Kronecker delta. In [2B] we will show how the general case corresponds to 
a "particle" moving in a non- commutative gauge field coupled with geodesic 
motion relative to the Levi-Civita connection associated with the g^. This 
result can be used to place the work of Tanimura [12] m a discrete context. 

It should be emphasized that all physics that we derive in this way is 
formulated in a context of non-commutative operators and variables. We 
do not derive electromagnetism, but rather a non-commutative analog. It 
is not yet clear just what these non- commutative physical theories really 
mean. Our initial idealisation of measurement is not the only model for 
measurement that corresponds to actual observations. Certainly the idea 
that we can measure time in a way that has "steps between the steps of 
time" is an idealisation. It happens to be an idealisation that fits a model of 
the universe as a cellular automaton. In a cellular automaton an observation 
is what an operator of the automaton might be able to do. It is not necessarily 
what the "inhabitants" of the automaton can perform. Here is the crux of the 
matter. The inhabitants can have only limited observations of the running of 
the automaton, due to the fact that they themselves are processes running on 
the automaton. The theories we build on the basis of DOC can be theories 
about the structure of these automata. They will eventually lead to theories 
of what can be observed by the processes that run on such automata. It 
is possible that the well known phenomena of quantum mechanics will arise 
naturally in such a context. These points of view should be compared with 

m. 

2.1 Brownian Walks and the Diffusion Equation 

To return to basics, consider the commutator equation in one space variable 
X. 
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[X, DX] = Jk 

for a single variable X. Written out, this equation becomes 

Jk = [X, J(X'-X)} = XJ(X'-X)-J(X'-X)X = J(X'(X'-X)-(X'-X)X). 

If k and the elements of the time series {X, X' , X" , ...} are all commuting 
scalars then this equation reduces to 

k= (X'-X) 2 . 

Thus 

X' = X± k 1/2 , 

a Brownian random walk, is a solution to the simplest one- dimensional com- 
mutator equation. 

Now lets examine this Brownian walk more closely by quantifying the 
time step as well as the space step. We take 

At = r 

so that 

DX = J(X' - X)/t 

where it is assumed that r is a scalar, commuting with all elements of the 
time series and commuting with the operator J (that is, r does not change 
with time). Now examine once again the equation 

[X, DX] = Jk. 

Let \X' — X\ — A. Then, repeating the calculation, we find 

k = (X' - Xf/r = A 2 jr. 

Hence 

A 2 /r = k. 

This tells us that if k is to be constant then there must be a constant rela- 
tionship between the square of the space interval for the Brownian walk and 
the size of the time interval. The remarkable point here is that it is just this 
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constant relationship that is required for a Brownian process to be described 
by the diffusion equation 

dP(x,t)/dt = Cd 2 P(x,t)/dx 2 

where the diffusion constant C is given by the formula 

C = A 2 /2r = k/2. 

The diffusion constant comes directly from our consideration involving the 
DOC commutator without any of the usual conceptual apparatus about ap- 
proximating a differential equation. 

To make this comparison, lets recall how the diffusion equation usually 
arises in discussing Brownian motion. We are given a Brownian process 
where 

x(t + r) = x(t) ± A 

so that the time step is r and the space step is of absolute value A. We regard 
the probability of left or right steps as equal, so that if P(x, t) denotes the 
probability that the Brownian particle is at point x at time t then 

P{x, t + r) = P(x- A, t)/2 + P(x + A)/2. 

From this equation for the probability we can write a difference equation for 
the partial derivative of the probability with respect to time: 

[{P(x, t + r)- P{x, t))/r] = (h 2 /2r) [{P{x - A, t) - 2P(x, t) + P(x + A))/A 2 ] 

The expression in brackets on the right hand side is a discrete approximation 
to the second partial of P(x, t) with respect to x. Thus if the ratio C = A 2 /2t 
remains constant as the space and time intervals approach zero, then this 
equation goes in the limit to the diffusion equation 

dP(x,t)/dt = Cd 2 P(x,t)/dx 2 . 

It is most curious how the diffusion constant comes up in these two con- 
texts. Lets try to think about the comparison between the non-commutative 
observational starting point and the more standard differential approxima- 
tion. In the non-commutative context we get A 2 from the appearance of the 
square of the difference of X' and X in the calculation of the commutator 
of X and DX. In the differential approximation, we get the A 2 from the ap- 
proximation of the second derivative of the probability P(x, t) with respect 
to x. The concept of probability does not appear in the non-commutative 
context. Clearly this subject needs more thought. 
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2.2 Planck's Numbers, Schrodinger's Equation and the 
Diffusion Equation 

First recall the Planck Numbers, h is Planck's constant divided by 2n. c is 
the speed of light. G is Newton's gravitational constant. The Planck length 
will be denoted by L, the Planck time by T and the Planck mass by M. Their 
formulas are 

M= ^hc/G 
L = h/Mc 
T = k/Mc 2 . 

These amounts of mass, length and time have just these dimensions and are 
constructed from the values of fundamental physical constants. They have 
roles in physics that point to deeper reasons than the formal for introducing 
them. Here we shall see how they are related to the Schrodinger equation. 

Recall that Schrodinger's equation can be regarded as the diffusion equa- 
tion with an imaginary diffusion constant. Recall how this works. The 
Schrodinger equation is 

ihdtp/dt = 

where the Hamiltonian H is given by the equation H = p 2 /2m + V where 
V(x,t) is the potential energy and p = h/id/dx is the momentum operator. 
With this we have p 2 /2m = (—h 2 /2m)d 2 /dx 2 . Thus with V(x,t) = 0, the 
equation becomes iUdtp/dt = {—h 2 /2m)d 2 il)/dx 2 which simplifies to 

d^/dt = {iTi/2m)d 2 %l)/dx 2 . 

Thus we have arrived at the form of the diffusion equation with an imagi- 
nary constant, and it is possible to make the identification with the diffusion 
equation by setting 

h/m = A 2 /t 

where A denotes a space interval, and r denotes a time interval as explained 
in the last section about the Brownian walk. With this we can ask what 
space interval and time interval will satisfy this relationship with a mass and 
Planck's constant? Remarkably, the answer is that this equation is satisfied 
when m is the Planck mass, A is the Planck length and t is the Planck time!! 
For note that 

L 2 /T = {h/Mc) 2 /(h/Mc 2 ) = h/M. 
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I now quote an email comment of Pierre Noyes: "With regard to your 
DOC derivation of the diffusion equation, and with an imaginary diffusion 
coefficient, the Schrodinger equation, note that the relation h/m = L 2 /T 
is satisfied for any mass m provided we take L = Compton wavelength = 
h/mc and T =Compton time= h/mc 2 — which is simply the time of a step 
length of this length taken at the velocity of light. I have a vague idea that I 
heard of this relation when I was a graduate student. In any case I am sure 
Feynman had it in mind when he used a random walk on the light cone to 
derive the 1 + 1 Dirac equation, and counted steps using z! So, in a sense, your 
DOC derivation of the diffusion equation does connect the Maxwell equations 
derivation via DOC, to the Dirac equation derivation — which in a vague 
sense was what I hoped we would be able to do this spring (2002). Of course 
this general result applies in particular to the Planck mass, which was your 
first observation. It is intriguing that if the mass scale is m [the Planck mass], 
then we can use either the Compton wavelength or the Schwarzchild radius 
at that mass scale as the step length in DOC. This reinforces my conviction 
(expressed long ago) that elementary particles are small black holes." 

The last part of Noyes' remark about the Schwarzchild radius refers to 
our work [37| explaining Ed Jones' microcosmology. Jones observed that if, 
for a particle of mass m we set the Schwarzchild radius (Rs = 2mG/c 2 ) equal 
to the Compton radius (Rc = h/2mc), then the resulting mass m is equal to 
one half the Planck mass! 

Rs — Rc 

2mG/c 2 = h/2mc 

m = (l/2)yJhc/G = M/2 

This is highly suggestive of limiting conditions on matter ( "Plancktonic mat- 
ter") prior to the Big Bang and leads in this way to specific cosmological 
predictions. It also gives an intriguing physical meaning to the Planck mass. 

What does all this say about the nature of the Schrodinger equation it- 
self? Interpreting it as a diffusion equation with imaginary constant suggests 
comparing with the DOC equation 

[X, DX] = JiC 

for a real constant C. This equation implicates a Brownian process where 
X' = X ± Z where Z 2 /t = iC. We can take Z = \fiL where L is a real step- 
length. This gives a Brownian walk in the complex plane with the correct 



15 



DOC diffusion constant. However, the relationship of this walk with the 
Schrodinger equation is less clear because the ip in that equation is not the 
probability for the Brownian process. To see a closer relationship we will 
take a different tack. 

Consider a discrete function if)(x,t) defined (recursively) by the following 
equation 

ip(x, t + r) = (i/2)ip(x - A, t) + (1 - i)i>(x, t) + (i/2)ip(x + A, t) 

In other words, we are thinking here of a random "quantum walk" where the 
amplitude for stepping right or stepping left is proportional to i while the 
amplitude for not moving at all is proportional to (1 — i). It is then easy to 
see that ip is a discretization of 

d^/dt = (iA 2 /2r)d 2 ^/dx 2 . 

Just note that ip satisfies the difference equation 

(i[>(x, t + r) -il>(x, t))/r = (iA 2 /2r)(^(x - A, f) - 2^{x, t) + if>(x + A, t))/A 2 

This gives a direct interpretation of the solution to the Schrodinger equa- 
tion as a limit of a sum over generalized Brownian paths with complex 
amplitudes. We can then reinterpret this in DOC terms by the equation 
[X, DX] = J(A 2 /r) or [X, DX] = 0, each of these contingincies happening 
probabilistically. It remains to be seen whether there is further insight to 
be gained into the Schrodinger equation via this combination of the DOC 
approach and the stochastic approach. 

2.3 DOC Chaos 

Along with the simple Brownian motion solution to the one dimensional com- 
mutator equation, there is a heirarchy of time series that solve this equation, 
with periodic and chaotic behaviour. These solutions can be obtained by 
taking 

X = J n Y 

where Y is a numerical scalar, and taking the commutator equation to be 

[X, DX\ = J 2n+1 k 
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where A; is a scalar. Expanding this equation, we find 

XJ(X' - X) - J(X' - X)X = J 2n+1 k 
J n YJ(J n Y' - J n Y) - J{J n Y' - J n Y)J n Y = J 2n+l k 

j2n+lyn+l /y> _ y^j — J 2 ™+ 1 ^y™+ 1 _ y n )Y = J 2n+1 k 

Y n+i(y> — Y) — {Y n+l - Y n )Y = k 

Y n + l (Y' -2Y) = k-Y n Y 

Y n+1 = (k-Y n Y)/(Y' -2Y). 

This last equation expresses the time series recursively where Y refers to the 
value of the series that is n time steps back from Y n . The first case of this 
recursion is 

Y" = (k-Y'Y)/(Y' -2Y). 

Next case is 

Y'" = {k-Y"Y)/(Y' -2Y). 

These recursions depend critically on the value of the parameter k. In the first 
case one sees periodic oscillations that (for appropriate values of k) destabilize 
and blow up, alternating between an unbounded phase and a bounded semi- 
periodic phase. We will investigate these time series in a separate paper. 



2.4 More Variables 

In the Feynman-Dyson derivation of electromagnetic formalism from com- 
mutation relations 24 one uses the relations 



[X t ,X 3 ] = 

[X i ,DX j ]=k8 ij 
where k is a scalar. Here we shall use 

[Xi,Xj] = 

[X u DXj] = JkSij 

as we did in analyzing the one-dimensional case. This allows us to have 
scalar evolution of the time series, but changes some of the issues in the 
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Feynman-Dyson derivation. These are in fact handled by the more general 
formalism that we discuss in the next two sections. Thus we shall aim in this 
section to see to what extent one can make simple models for this version 
of the Feynman-Dyson relations. Models of this sort will be another level of 
approximation to discrete electromagnetism. 

Writing out the commutation relation [X, DX] = Jk, and not making 
any assumption that X' commutes with X, we find 

J~ l [X, DX] = X\X' - X) - (X' - X)X 

= X\X' - X) - X(X' -X)+ X(X' - X) - (X' - X)X 

= (X' - X) 2 + (XX' - X'X) = (X' - X) 2 + [X, X'}. 

Thus the commutation relation [X, DX] = Jk becomes the equation 

(X'-X) 2 + [X,X'\ =k. 

By a similar calculation, the equation [X, DY] = becomes the equation 

(X' - X)(Y' -Y) + [X, Y'] = 0. 

These equations are impossible to satisfy simultaneously for k ^ if we 
assume that X and X' commute and that X and Y' commute and that 
[Y, DY] = Jk. For then we would need to solve: 

(X' - X) 2 = k. 

{Y 1 - Y) 2 = k. 

(X 1 - X)(Y' -Y) = 0. 

with the first two equations implying that (X — X') and (Y — Y') are each 
non-zero, and the third implying that their product is equal to zero. In other 
words, the equations below cannot be satisfied if the time series are composed 
of commuting scalars. 

[X, DX] = Jk 
[Y, DY] = Jk 
[X,Y]=0 

In order to make such models we shall have to introduce non-commutativity 
into the time series themselves. In a certain sense this is analoguous to 
the introduction of non-commutative algebra in the Dirac equation in 3 + 
1 dimensions, and to the introduction of non-commutative fields in gauge 
theory. 
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Here is an example of such a model. 
Return to the equations 

(X'-X) 2 + [X,X'\ = k. 

(X' - X){Y' -Y) + [X, Y'} = 

expressing the behaviour for two distinct variables X and Y. If [X, X'] = 0, 
then we have (X' — X) 2 = k so that 

X' = X ± Vk. 

In order for the second equation to be satisfied, we need that 

[X, Y'\ = ±k 

where the ambiguity of sign is linked with the varying signs in the temporal 
behaviour of X and Y. We will make the sign more precise in a moment, but 
the radical part of this suggestion is that for two distinct spatial variables X 
and Y, there will be a commutation relation between one and a time shift of 
the other. 

If the space variables are labeled Xj, then we can write 

X* +1 = X\ + e\k 

where e" is plus one or minus one. Thus each space variable performs a walk 
with the fixed step-length k. We shall write informally 

X[ = X + tik 

where it is understood that the epsilon without the superscript connotes the 
sign change that occurs in this juncture of the process. We then demand the 
commutation relations 

[X' i ,X j ] = [X' j ,X i ] = ei e j k. 

Each Xj is a scalar in its own domain, but does not commute with the time 
shifts of the other directions. We then can have the full set of commutation 
relations: 

[X' i ,X j ] = [X' j ,X i ] = e i e j k. 
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[X t ,X 3 } = 
[X i ,DX j ] = Jk5 ij 

so that the system will satisfy the assumptions supporting the Feynman- 
Dyson derivation. In this system,the elements of a given time series Xj, X-, X", ■ ■ ■ 
commute with one another. The basic field element in the Feynman-Dyson 
set up is the magnetic field B defined by the (non-commutative) vector cross 
product 

B = (l/k)DX x DX. 

Here we have 

DXi = J{X[ - Xi) = JeiVk. 

Thus 

B = J 2 e' x e 

where e = (ej, e 2 , e 3 ) (assuming three spatial coordinates) and e' denotes this 
vector of signs at the next time step. In this way we see that we can think 
of each spatial coordinate as providing a long temporal bit string and the 
three coordinates together give the field in terms of the vector cross product 
of their temporal cross sections at neighboring instants. It is interesting to 
compare this model with the color algebra in the following paper by Wene 

03- 

2.5 Discrete Classical Electromagnetism 

It is of interest to compare these results with a direct discretization of classical 
electromagnetism. Suppose that X, X', X", X'", • ■ ■ is a time series of vectors 
in R 3 (where R denotes the real numbers). Let dX = X' — X be the usual 
discrete derivative (with time step equal to one for convenience). Let A* B 
denote the usual inner product of vectors in three dimesions. Assume that 
there are fields E and B such that 

d 2 X = E + dX x B 

(the Lorentz force law). Assume also that E and B are perpendicular to the 
velocity vector dX, and that E is perpendicular to B. Then we have 

dX' x dX = (dX' - dX) x dX = (d 2 X) x (dX) 
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= E x dX + (dX x B) x dX 

= ExdX - dX(B • dX) + (dX • dX)B. 

Since E is perpendicular to dX we know there is a A such that E x dX = \B 
and we have B • dX = since B is perpendicular to dX. Therefore 

dX' x dX = XB + \\dX\\ 2 B 

so that 

B = dX' x rfX/(A+ ||rfX|| 2 ). 

Up to the factor in the denominator, this formula is in exactly the same 
pattern as the formula in our discrete model for DOC electromagnetism as 
described in the previous subsection. To see this, note that the B field in 
the DOC model is proportional to DX x DX and that DX = JdX so that 
DX x DX = JdX x JdX = J 2 dX' x dX. Up to the time-shifting algebra and 
a proportionality constant, the expressions are the same! Clearly more work 
is needed in comparing classical discrete electromagnetism with the results 
of a discrete analysis of the Feynman-Dyson derivation. 

3 Gauge Fields and Differential Geometry 

Letting X i [i = 1, 2, d) denote a set of spatial variables (non- commutative 
time series in the sense of our discrete ordered calculus), we will look at 
a collection of basic assumptions about the commutation of these variables 
and of their derivatives. It is natural from the point of view of the discrete 
ordered calculus to have 

[X, t ,X 3 ] = 

for all i and j. There are no other natural commutations from the point of 
view of this calculus. 

We shall define by the equation 

[X h Xj] = gij. 
Here Xj is shorthand iorDXj and 

[A, B] — AB — BA. 
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Along with this commutator equation, we will assume that 



[Xi,Xj] 







and 




0. 



Here it is assumed that gij is non-degenerate in the sense that there exists 
g tj so that 



Here we are using the Einstein summation convention that implicitly assumes 
that we sum over repeated indices in an expression. Symbol <5] is a Kronecker 
delta, equal to 1 when i equals j and otherwise. 

The first result that is a direct consequence of these assumptions is the 
symmetry of the "metric" coefficients g % K That is, we shall show that 



9 J 9jk 



and that 



9ij9 



5l 



Lemma 3. g {j = g^. 



Proof. 



9ij 9ji 



[Xi,Xj] - [Xj,Xi] 
[Xi, Xj] + [Xi, Xj] 
— D[Xi, Xj\ 



0. 



For the purpose of doing calculus in this situation we define X i by the 
equation 
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X* = g lk X k . 

The operator X % is simply the index shift of the corresponding Xj. We do 
not define a corresponding X 1 . It is easy to check the equation 

Consequently we define the derivative of an operator F with respect to Xi 
by the equation 

9*F = [F, X*\ 

and the corresponding lowered derivative by the formula 

diF=[F,Xi]. 

Note that we have 

diXj = Qij. 

We also define 

diF=[Xi,F], 

the derivative of F with respect to the conjugate variable X 1 . 
With these partial derivatives in hand, we define F by the formula 

F = d k FX k . 

If F commutes with g' lj then it is easy to see that 

F = d k FX k . 

These formulas extend (implicitly) the definition of the time series to entities 
other than the operators Xj since 

p = DF = J(F' - F). 

A stream of consequences then follows by differentiating both sides of the 
equation 
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Note that 



9ij — [Xi,Xj\. 



g\j = [X i ,X j ] + [X i ,D 2 X j ] 

by the Leibniz rule 

D[A,B] = [DA.B] + [A,DB]. 
Note also that we can freely use the Jacobi identity 

[A, [B,C]] + [C, [A,B]] + [B,[C,A]] = 0. 

In particular, the Levi-Civita connection 

T ijk = (l/2)(dig jk + djg ik - d k g i:j ) 

associated with the g i3 - comes up almost at once from the differentiation pro- 
cess described above. To see how this happens, view the following calculation 
where 

d l d 3 F=[X l ,[X 3 ,F]]. 

We apply the operator fydj to the second DOC derivative of X k . 

Lemma 4. T ijk = (l/2)d i d j D 2 X k 
Proof. 

d t d J D 2 X k = [X t ,[X„D 2 X k ]] 
— [Xi,gj k — [Xj,X k ]\ 
= [Xi,g'jk\ - [Xi, [Xj,X k ]] 
— [Xi,gjk] + [Xk, [Xi,Xj]\ + [Xj, [X k ,Xi] 
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— [djk, Xi] + [X k , gij] + [Xj, —gik] 



— diQjk — dk9ij + djgik 

It is remarkable that the form of the Levi-Civita connection comes up directly 
from this non-commutative calculus without any apriori geometric interpre- 
tation. We shall discuss the context of this result in the next two sections of 
the paper. 

One finds that 

D 2 X t (!, ■ //,,,(/,. F'\\-.' + T tjk XJX k 

where 

F rs = [X r ,X s ]. 
It follows from the Jacobi identity that 

satisfies the equation 

diFjk + 9jF ki + dtFij = 0, 

identifying F^ as a non- commutative analog of a gauge field. Gj is a non- 
commutative analog of a scalar field. The details of these calculations will 
be found in 

This description of the equations for a non- commutative particle in a 
metric field illustrates the role of the background discrete time in this theory. 
In terms of the background time the metric coefficients are not constant. It 
is through this variation that the spacetime derivatives of the theory are 
articulated. The background is a process with its own form of discrete time, 
but no spacetime structure as we know and observe it. Our observation of 
spacetime structure appears as a rough (commutative) approximation to the 
processes described as consequences of the basic non-commutative equations 
of the discrete ordered calculus. 
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4 Curvature, Jacobi Identity and the Levi- 
Civita Connection 



In this section, we go back to basics and examine the context of calculus 
defined via commutators. We shall use a partially index-free notation. In 
this notation, we avoid nested subscripts by using different variable names 
and then using these names as subscripts to refer to the relevant variables. 
Thus we write X and Y instead of Xi and Xj, and we write gxy instead 
of gij. It is assumed that the derivation DX has the form DX = [X, J] for 
some J. 

The bracket [A, B] is not assumed to be a commutator. It is assumed to 
satisfy the Jacobi identity, bilinearity in each variable, and the Leibniz rule 
for all functions of the form 5k (A) = [A, K\. That is we assume that 



Recall that in classical differential geometry one has the notion of a co- 
variant dervative, defined by taking a difference quotient using parallel trans- 
lation via a connection. Covariant derivatives in different directions do not 
necessarily commute. The commutator of covariant derivatives gives rise to 
the curvature tensor in the form 



If derivatives do not commute then we regard their commutator as express- 
ing a curvature. In our non-commutative context this means that curvature 
arises 'prior to any notion of covariant derivatives since even the basic deriva- 
tives do not commute. 

We shall consider derivatives in the form 



S K (AB) = 5 K (A)B + A5 K (B). 



[V < ,V J ■]X fc = i^J J .X , . 



V X (A) 



[A, A x ]. 



Examine the following computation 



VxVyF= [[F,A Y ],A X ] 



[[A x ,F},A Y }-[[A Y ,Ax},F] 



[[F,A X },A Y } + [[A X ,A Y },F] 
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= V Y V X F+[[A X ,A Y ],F]. 

Thus 

[V X ,V Y ]F = R XY F 

where 

RxyF = [[A x ,A y ],F]. 
We can regard R XY as a curvature operator. 

The analog in this context of flat space is abstract quantum mechanics! 
That is, we assume position variables (operators) X, Y, ■ ■ ■ and momentum 
variables (operators) P X) P Y) ■ ■ ■ satisfying the equations below. 

[X,Y]=0 

[P x ,P Y ] = 
[X,P Y }=5 XY 

where 5 XY is equal to one if X equals Y and is zero otherwise. We define 

d x F = [F,P X ] 

and 

d Px F = [X,F}. 

In the context of the above commutation relations, note that these derivatives 
behave correctly in that 

d x (Y) =5 XY 

and 

d Px (P Y ) = 5 XY 

d Px (Y)=0 = d x (P Y ) 

with the last equations valid even if X — Y. Note also that iterated par- 
tial derivatives such as d x d Y commute. Hence the curvature R XY is equal 
to zero. We shall regard these position and momentum operators and the 
corresponding partial derivatives as an abstract algebraic substitute for flat 
space. 
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With this reference point of (algebraic, quantum) flat space we can define 

P x = Px- Ax 

for an arbitrary algebra-valued function of the variable X. In indices this 
would read 

Pi Pi A{, 

and with respect to this deformed momentum we have the covariant deriva- 
tive 

V X F = [F, Py] = [F, Py + Ay] = dyF + [F, Ay}. 

The curvature for this covariant derivative is given by the formula 

RxyF = [Vx, Vy]F = [[A X , Ay], F] 
where \x = Px — Ax- Hence 

Rxy = [Px - A x , P Y - A Y ] = -[P x , Ay] - [A x , Py] + [A x , Ay] 

= d X Ay - d Y A X + [A X , Ay}. 

With indices this reads 

Rij = d{Aj djA,i + [Ai, Aj]. 

and the reader will note that this has the abstract form of the curvature of a 
Yang-Mills gauge field, and specifically the form of the electromagnetic field 
when the potentials Ai and Aj commute with one another. 

Continuing with this example, we compute 

[X, Py] = [X, Py - Ay] = 5 XY - [X , Ay}. 

Let 

9xy = $xy - [X, Ay] 

so that 

[X,Py]=g X Y- 

We will shortly consider the form of this general case, but first it is useful to 
restrict to the case where [X, Ay] = so that qxy = <$xy- This is the domain 
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to which the original Feynman-Dyson derivation applies. In order to enter 
this domain, we set 

X = DX = P X = P X - A x . 

We then have 

[Xi,Xj] = 
[Xi,Xj] = Sij 

and 

Rij = [Xi, Xj] = diAj - djAi + [Ai, Aj\. 

Note that even under these restrictions we are still looking at the possibility 
of a non-abelian gauge field. The pure electromagnetic case is when the 
commutator of Ai and Aj vanishes. But why do we set X = P x l The 
answer to this is the key to the gauge interpretation of electromagnetism, 
for with this interpretation we find that X satisfies the Lorentz force law 
X = E + X x B where B represents the magnetic field and E the electric 
field (in the case of three space variables Xi with i — 1,2, 3.) To see how this 
works, suppose that Xi = Ei + EijXj and suppose that Ei and commute 
with X k . Then we can compute 

[Xi, Xj] = [X, h Ej + Fj k X k ] 

= Fj k [Xi,X k ] = Fj k 5 ik = Fji. 

This implies that 

Fij = [X h Xj] = R l3 = diAj - djAi + [A i} Aj] 

since [Xj,Xj] + [Xj,Xj] = D[Xi,Xj] = 0. It is then easy to verify that the 
Lorentz force equation is satisfied with B k = €ij k Rij and that in the case of 
[Ai, Aj] = this leads directly to standard electromagnetic theory when the 
bracket is a Poisson bracket (see the next section for a discussion of Poisson 
brackets). When this bracket is not zero but the potentials Ai are functions 
only of the Xj we can look at a generalization of gauge theory where the non- 
commutativity comes from internal Lie algebra parameters. This shows that 
the Feynman-Dyson derivation supports certain generalizations of classical 
electromagnetism, and this will be the subject of a more expanded version 
of this paper. 
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In regard to this last remark, the reader should note that in our |28| 
127] algebraic and discrete version of the Feynman-Dyson derivation it was 
actually an additional assumption that B x B = where B x B denotes the 
(non- commutative) vector cross product of B with itself. (Note that B = 
(1/2)X xX .) In the original Dyson paper this cross product vanished because 
of assumptions about the operators and their Hilbert space representations. 
With B x B as an extra term, the Feynman-Dyson derivation is indeed a 
non-commutative generalization of electromagnetism and includes forms of 
gauge theories among its models. 

Generalizing, we wish to examine the structure of the following special 
axioms for a bracket. 

[X,DY] =g XY 
[X,Y] = 
[Z,g XY }=0 
[gxY,gzw] = o 

Note that 

Dg YZ = D[Y, DZ] = [DY, DZ] + [Y, D 2 Z]. 
and that D [X, g X y] = implies that 

[g XY ,DZ] = [Z,Dg XY }. 
Define two types of derivations as follows 

V X (F) = [F,DX] 

and 

W DX (F) = [X,F]. 

These are dual with respect to g XY and will act like partials with respect to 
these variables in the special case when g XY is a Kronecker delta, 5 XY . If the 
form g XY is invertible, then we can rewrite these derivations by contracting 
the inverse of g to obtain standard formal partials. 
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V DX V DY D 2 Z = [X, [Y,D 2 Z]] 
= [X, Dg YZ - [DY, DZ]] = [X, Dg YZ ] - [X, [DY, DZ\\ 
= [g YZ ,DX]-[X, [DY, DZ}} 
= V x (gyz)-[X, [DY, DZ]]. 
Now use the Jacobi identity on the second term and obtain 

V DX V DY D 2 Z = V x (g YZ ) + [DZ, [X,DY]] + [DY, [DZ,X\] 

= Vx(fl'yz) - Vz(s-xy) + Vy(#xz)- 
This is the formal Levi-Civita connection. 

At this stage we face once again the mystery of the appearance of the 
Levi-Civita connection. There is a way to see that the appearance of this 
connection is not an accident, but rather quite natural. We shall explain 
this point of view in the next section where we discuss Poisson brackets 
and the connection of this formalism with classical physics. On the other 
hand, we have seen in this section that it is quite natural for curvature in 
the form of the non-commutativity of derivations to appear at the outset 
in a non-commutative formalism. We have also see that this curvature and 
connection can be understood as a measurement of the deviation of the the- 
ory from the "flat" commutation relations of ordinary quantum mechanics. 
Electromagnetism and Yang-Mills theory can be seen as the theory of the 
curvature introduced by such a deviation. On the other hand, from the 
point of view of metric differential geometry, the Levi-Civita connection is 
the unique connection that preserves the inner product defined by the metric 
under the parallel translation defined by the connection. We would like to 
see that the formal Levi-Civita connection produced here has this property 
as well. 

To this end lets recall the formalism of parallel translation. The infinites- 
imal parallel translate of A is denoted by A' = A + 5 A where 

5A k = -T^AUX 1 

where here we are writing in the usual language of vectors and differentials 
with the Einstein summation convention for repeated indices. We assume 
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that the Christoffel symbols satisfy the symmetry condition T k - = T^. The 
inner product is given by the formula 

< A, B >= 9ij A l Bi 

Note that here the bare symbols denote vectors whose coordinates may be 
indicated by indices. The requirement that this inner product be invariant 
under parallel displacement is the requirement that 5{gijA l A>) = 0. Calcu- 
lating, one finds 

SigijA'A*) = (d k g tj )A l A^dX k + g i3 5(A l )A^ + gij A%Ai) 

= {d k g ij )A i A^dX k - 9ij r rs A r dX s Ai - g l3 A^ rs A T dX* 
= {d k g ij )A i AUX k - ;/,,r;,. r.FV/A- - g,^\.Y.VdX* 
= (d k g tj )A*AidX k - g sj T s ik A i A j dX k - gJ^.Y A'd.X 1 

Hence 

(d k gij) = g sj r s ik + g is r s jk . 

From this it follows that 

r ijk = gis^jk = (l/2)(d k g ij - dig jk + d 3 (g ik )). 

Certainly these notions of variation can be imported into our abstract con- 
text. The question remains how to interpret the new connection that arises. 
We now have a new covariant derivative in the form 

V t X j = diX j + T j kl X k . 

The question is how the curvature of this connection interfaces with the 
gauge potentials that gave rise to the metric in the first place. The theme 
of this investigation has the flavor of gravity theories with a qauge theoretic 
background. We will investigate these relationships in detail in a sequel to 
this paper. 



32 



5 Poisson Brackets and Commutator Brack- 
ets 

Dirac JT] introduced a fundamental relationship between quantum mechan- 
ics and classical mechanics that is summarized by the maxim replace Poisson 
brackets by commutator brackets. Recall that the Poisson bracket {A, B} is 
defined by the formula 

{A, B} = (dA/dq)(dB/dp) - (dA/dp)(dB/dq), 
where q and p denote classical position and momentum variables respectively. 

In our version of discrete physics the noncommuting variables are func- 
tions of discrete time, with a DOC derivative D as described in the first 
section. Since DX = X J — JX = [X, J] is itself a commutator, it follows 
that 

D{[A,B)) = [DA, B] + [A,DB] 

for any expressions A, B in our ring R. A corresponding Leibniz rule for 
Poisson brackets would read 

(d/dt){A, B} = {dA/dt, B} + {A, dB/dt}. 

However, here there is an easily verified exact formula: 

(d/dt){A, B} = {dA/dt, B} + {A, dB/dt} - {A, B}(dq/dq + dp /dp). 

This means that the Leibniz formula will hold for the Poisson bracket exactly 
when 

(dq/dq + dp /dp) = 0. 

This is an integrability condition that will be satisfied if p and q satisfy 
Hamilton's equations 

q = dH/ dp, 
p = —dH/dq. 
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This, of course, means that q and p are following a principle of least ac- 
tion with respect to the Hamiltonian H . Thus we can interpret the fact 
D([A,B]) = [DA, B] + [A, DB\ in the discrete (commutator) context as an 
analog of the principle of least action. Taking the discrete context as funda- 
mental, we say that Hamilton's equations are motivated by the presence of 
the Leibniz rule for the discrete derivative of a commutator. The classical 
laws are obtained by following Dirac's maxim in the opposite direction! Clas- 
sical physics is produced by following the correspondence principle upwards 
from the discrete. 

Taking the last paragraph seriously, we must reevaluate the meaning of 
Dirac's maxim. The meaning of quantization has long been a basic mystery 
of quantum mechanics. By traversing this territory in reverse, starting from 
the non- commutative world, we begin these questions anew. 

In making this backwards journey to classical physics we see how our ear- 
lier assertion that bare quantum mechanics of commutators can be regarded 
as the background for the coupling with other fields (as in the description of 
formal gauge theory in the last section), fits with Poisson brackets. The bare 
Poisson brackets satisfy 

UhQj} = 
{PhPj} = o 

{<lhpj} = Sij- 

In our previous formalism, we would identify Xj as the correspondent with qi 
and Pj as the correspondent of pj. And, given a classical vector potential A, 
we could write the coupling dqi/dt = pi — Ai to describe the motion of a par- 
ticle in the presence of an electromagnetic field. The analog of the Feynman 
Dyson derivation is then expressed classically in terms of the Poisson brack- 
ets. Similar remarks apply to the analogs for gauge theory and curvature. 
In particular it is of interest to see that our derivation of the Levi-Civita 
connection corresponds to the motion of a particle in generalized coordinates 
that satisfies Hamilton's equations. The fact that such a particle moves in a 
geodesic according to the Levi-Civita connection is a classical fact that was 
surely one of the motivations for the development of differential geometry. 
Our derivation of the Levi-Civita connection, interpreted in Poisson brackets, 
reproduces this result. 
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To see how this works, let ds 2 = g lJ dxjdxj denote the metric in the 
generalized coordinates x^. Then the velocity of the particle has square v 2 = 
(ds/dt) 2 = g^XiXj. The Lagrangian for the system is the kinetic energy 
L = mv 2 /2 = mg^XiXj/2. Then the canonical momentum is pj = dL/dxj, 
and with qi = xi we have the Poisson brakets 

= UhPj} = {xudL/dxj} = {x i ,mg 3k x k }. 

Taking m = 1 for simplicity, we can rewrite this bracket as 

This, in Poisson brackets, is our generalized equation of motion. 

The classical derivation applies Lagrange's equation of motion to the sys- 
tem. Lagrange's equation reads 

djdtipLjdxi) = dL/dxi. 

Since this equation is equivalent to Hamilton's equation of motion, it follows 
that the Poisson brackets satisfy the Leibniz rule. With this, we can proceed 
with our derivation of the Levi-Civita connection in relation to the accelera- 
tion of the particle. In the classical derivation, one writes out the Lagrange 
equation and solves for the acceleration. The advantage of using only the 
Poisson brackets is that it shows the relationship of the connection with the 
Jacobi identity and the Leibniz rule. 

This discussion raises further questions about the nature of the general- 
ization that we have made. Originally Hermann Weyl j33] generalized classi- 
cal differential geometry and discovered gauge theory by allowing changes of 
length as well as changes of angle to appear in the holonomy. Here we arrive 
at a very similar situation via the properties of a non- commutative discrete 
calculus of observations. A closer comparison with the geometry of gauge 
theories is called for. 



6 Discussion on ^-Deformation 

The direct relation between the content of local physical descriptions based on 
the DOC calculus and more global considerations are a matter of speculation. 
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One strong hint is contained in the properties of the discrete derivative that 
has the form 

D q f( x ) = (fM - f(x))/(qx - x). 
The classical derivative occurs in the limit as q approaches one. 

In the setting of q not equal to one, the derivative D q is directly related 
to fundamental noncommutativity. Consider variables x and y such that 
yx = qxy where q is a commuting scalar. Then the expansion of (x + y) n 
generates a g-binomial theorem with g-choice coefficients composed in q- 
factorials of g-integers [n] q where 

[n] q = l + q + q 2 + ... + q^- 1 \ 

The derivative D q is directly related to the g-integers via the formula 

D q {x n ) = [n} q x n -\ 



In the context of this paper, we have considered discrete derivatives in 
the form 

d A f(x) = (f(x + A)-f(x))/A. 
This will convert to the g-derivative if x + A = qx. Thus we need 

g = (x + A)/x. 

This means that a direct translation from DOC to g-derivations could be 
effected if we allowed g to vary as a function of x and introduced the temporal 
operator J into the calculus of g-derivatives. 

In general, many g-deformed structures such as the quantum groups as- 
sociated with the classical Lie algebras appear to be entwined with the dis- 
cretization inherent in D q . The quantum groups have turned out to be deeply 
connected with topological amplitudes for networks describing knots and 
three dimensional spaces. (See the next section of this paper.) The ana- 
log for the quantum groups in dimension four is being sought. If there is 
a connection between the local and the global parts of our essay it may lie 
in hidden connections between discretization and quantum groups. Clearly 
there is much work to be done in this field. 
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There is a clue about the meaning of the operator J (DF = [F, J] in 
the discrete ordered calculus) in the context of quantum groups. Quantum 
groups are Hopf algebras. A quantum group such as G = U q (SU(2)) is 
actually an algebra over a field k with an antipode 



S : G >G 



and a coproduct 

A : G — > G <g> G, 

a unit 1 and a couinit 

e : G — > k. 

The coproduct is a map of algebras. The antipode is an antimorphism, 
S(xy) = S(y)S(x), and generalizes the inverse in a group in the sense that 
T,S(xi)X2 = e(x)l and T l xiS(x 2 ) = e(x)l where A(x) = T,x± <g> x 2 . 

An element g in a quantum group G is said to be a grouplike element 
if A(g) = g ® g and S(g) = g~ l . In many quantum groups (such as G = 
C/ g (577(2))) the square of the antipode is represented via conjugation by a 
special grouplike element that we shall denote by J. Thus 

s 2 ( x ) = r l x.j 

for all x in G. This means that it is possible to define the discrete ordered 
calculus in the context of a quantum group G (as above) by taking J to be 
the special grouplike element. Then we have 

DX = [X, J}= XJ - JX = J(J- l XJ -X) = J(S 2 (X) - X). 

Conjugation by the special grouplike element in the quantum group consti- 
tutes the time evolution operator in this algebra. 

There are a number of curious aspects to this use of the discrete ordered 
calculus in a quantum group. First of all, it is the case that in some quantum 
groups (for example with undeformed classical Lie algebras) the square of the 
antipode is equal to the identity mapping. From the point of view of DOC, 
time does not exist in these algebras. But in the g-deformations such as 
U q (SU(2)), the square of the antipode is quite non-trivial and can serve well 
as the tick of the clock. In this way, g-deformations do provide a context for 
time. In particular, this suggests that the q- deformations of classical spin 
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networks [3H] should be able to accommodate time. A suggestion directly 
related to this remark occurs in jU], and we shall take this up at the end of 
the next section of this paper. 

7 Networks, Discrete Spacetime and the Dirac 
Equation 

One can consider replacing continuous space (such as Euclidean space with 
the usual topology) by a discrete structure of relationships. The geometry of 
the Greeks held a discrete web of relationships in the context of continuous 
space. That space was not coordinatized in our way, nor was it held as an 
infinite aggregate of points. In general topology there is a wide choice for 
possible spatial structures (where we mean by a space a topology on some 
set). 

Discretization of space and time implicates the replacement of spacetime 
by a network, graph or complex that has nodes for the points and edges to 
indicate significant relationships among the points. 

Euler's work in the eighteenth century brought forth the use of abstract 
graphs as holders of spatial structure. After Euler it was possible to find the 
classification of the Greek regular solids in the the (wider) classification of 
the regular graphs on the surface of the sphere. Metric can disappear into 
relationship under the topological constraint of Euler's formula V—E+F = 2, 
where V denotes the number of vertices, E the number of edges and F the 
number of faces for the connected graph G on the sphere. 

A network itself can represent an abstract space. Embeddings of that 
network into a given space (such as graphs on the two dimensional sphere) 
correspond to global constraints on the structure of the abstract graph. 

Now a new theme arises, motivated by a conjunction of combinatorics and 
physics. Imagine labelling the edges of the network from some set of "colors" . 
These colors can represent the basic states of a physical system, or they can 
be an abstract set of distinct markers for purely mathematical purposes. 
Once the network is labelled, each vertex is an entity with a collection of 
labels incident to it. Let there be given a function that associates a number 
(or algebra element) to each such labelled vertex. Call this number the vertex 
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weight at that vertex. Let C denote a specific coloring of the network iV and 
consider the product, over all the vertices of N of the values of the vertex 
weights. Finally let Z(N) , the amplitude of the network, be defined as the 
summation of the product of the vertex weights over all colorings of the net. 
Z(N) is also called the partition function of the network. 

Amplitudes of this sort are exactly what one computes in finding the par- 
tition function of a physical system or the quantum mechanical amplitude 
for a discrete process. In all these cases the network is interwoven with the 
algebraic structure of the vertex weights. It is only recently that topological 
properties of networks in three dimensional space have come to be under- 
stood in this way |22], [Tj . |45j . This has led to new information about the 
topology of low dimensional spaces, and new relationships between physics 
and topology. 

A classical example of such an amplitude was discovered by Roger Penrose 
jS] in elucidating special colorings of 3-regular graphs in the plane. A 3- 
regular graph G has three edges incident to each vertex. When embedded in 
the plane, these edges acquire a specific cyclic order. Three colors are used. 
One associates to each vertex the weight 



where a,b,c denote the edges meeting the vertex in this cyclic order, and 
the epsilon is equal to 1, —1 according as the edges have distinct labels in 
the given or reverse cyclic order, or if there is a repetition of labels. The 
resulting amplitude counts the number of ways to color the network with 
three colors so that three distinct colors are incident to each vertex. This 
result is a perspicuous generalization of the classical four color problem of 
coloring maps in the plane with four colors so that adjacent regions receive 
different colors. 

The Penrose example generalizes to networks whose amplitudes embody 
geometrical properties of Euclidean three dimensional space (angles and their 
dependence). Geometry begins to emerge in terms of the averages of prop- 
erties of an abstract and discrete network of relationships. Topological prop- 
erties emerge in the same way. The idea of space may change to the idea of 
a network with global states and a functor that associates this network and 
its states to the more familiar properties that a classical observer might see. 
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7.1 Remarks on Quantum Mechanics 

We should remark on the basic formalism for amplitudes in quantum me- 
chanics. The Dirac notation (A\B) [TI] denotes the probability amplitude 
for a transition from A to B. Here A and B could be points in space (for 
the path of a particle), fields (for quantum field theory), or geometries on 
spacetime (for quantum gravity). The probability amplitude is a complex 
number. The actual probability of an event is the absolute square of the 
amplitude. If a complete set of intermediate states Ci,C2,...C n is known, 
then the amplitude can be expanded to a summation 

(A\B) = Y!U{A\C i ){C i \B). 

This formula follows the formalism of the usual rules for probability, and it 
allows for the constructive and destructive interference of the amplitudes. It 
is the simplest case of a quantum network of the form 

A * C * B 

where the colors at A and B are fixed and we run through all choices of 
colors for for the middle edge. The vertex weights at the vertices labelled 
* are (A\C) and (C\B) respectively. A measurement at the C edge reduces 
the big summation to a single value. 

Consider the generalization of the previous example to the graph 

A * C 1 * C 2 * ... * C m B 

With A and B fixed the amplitude for the net is 

< A\B >= Vi&K^Kn < A\C\ >< Cl\C\ > ... < C?JB > 

One can think of this as the sum over all the possible paths from A to 
B. In fact in the case of a "particle" travelling between two points in space, 
this is exactly what must be done to compute an amplitude - integrate over 
all the paths between the two points with appropriate weightings. In the 
discrete case this sort of summation makes perfect sense. In the case of 
a continuum there is no known way to make rigorous mathematical sense 
out of all cases of such integrals. Nevertheless, the principles of quantum 
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mechanics must be held foremost for physical purposes and so such "path 
integrals" and their generalizations to quantum fields are in constant use by 
theoretical physicists ^H] who take the point of view that the proof of a 
technique is in the consistency of the results with the experiments. When 
the observations themselves are mathematical (such as finding invariants of 
knots and links), the issue acquires a new texture. 

Now consider the summation discussed above in the case where n = 2. 
That is, we shall assume that each C fc can take two values, call these values 
L and R. Furthermore let us suppose that < L\R >=< R\L >= ^f^l while 
< L\L >=< R\R >= 1. The amplitudes that one computes in this case 
correspond to solutions to the Dirac equation ^1] in one space variable and 
one time variable. This example is related to an observation of Richard 
Feynman JH]. in [25J we give a very elementary derivation of this result 
and we show how these amplitudes give solutions to the discretized Dirac 
equation, so everything is really quite exact and one can understand just what 
happens in taking the limit to the continuum. In this example a state of the 
network consists in a sequence of choices of L or R. These can be interpreted 
as choices to move left or right along the light-cone in a Minkowski plane. 
It is in summing over such paths in spacetime that the solution to the Dirac 
equation appears. In this case, time has been introduced into the net by 
interpreting the sequence of nodes in the network as a temporal direction. 

More specifically, let (a, b) denote a point in discrete Minkowski spacetime 
in lightcone coordinates. This means that a denotes the number of steps 
taken to the left and b denotes the number of steps taken to the right. We 
let ipL(a,b) denote the sum over the paths that enter the point (a, b) from 
the left and ipn{a, b) the sum over the paths that enter (a, b) from the right. 
Each path P contributes i c ^ p where c(P) denotes the number of corners in 
the path. View the diagram below. 



(a,b+l) 



(a,b) 
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It is clear from the diagram that 

i/; L (a, b + 1) = i> L (a, b) + ii/i R (a, b). 

Thus we have that 

di> L /dR = iip R 

and similarly 

dip R /dL = iijj L . 

This pair of equations is the Dirac equation in light cone coordinates. 

This discrete derivation of the Dirac equation is simpler than the method 
used in |25j . I am indebted to Charles Bloom [3| for pointing this out to 
me. In fact, this form of the discretization is essentially Feynman's original 
method as is evident from the reproduction of Feynman's handwritten notes 
in Figure 8 of the review paper by Schweber. For one approach, very 
close in spirit, that generalizes this exercise of Feynman to four dimensional 
discrete spacetime see |4T)] . 

As in the Dirac equation example, one way to incorporate spacetime is 
to introduce a temporal direction into the net. At a vertex, one must specify 
labels of before and after to each edge of the net that is incident to that 
vertex. If there is a sufficiently coherent assignment of such local times, then 
a global time direction can emerge for the entire network. Networks endowed 
with temporal directions have the structure of morphisms in a category where 
each morphism points from past to future. A category of quantum networks 
emerges equipped with a functor (via the algebra of the vertex weights) 
to morphisms of vector spaces and representations of generalized symmetry 
groups. Appropriate traces of these morphisms produce the amplitudes. 

Quantum non-locality is built into the network picture. Any observer 
taking a measurement in the net has an effect on the global set of states 
available for summation and hence affects the possibilities of observations at 
all other nodes in the network. By replacing space with a network we obtain 
a precursor to spacetime in which quantum mechanics is built into the initial 
structure. 

Remark. A striking parallel to the views expressed in this section can be 
found in [T2]. Concepts of time and category are discussed by Louis Crane 
jH], [SI in relation to topological quantum field theory. In the case of Crane's 
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work there is a deeper connection with the methods of this paper, as I shall 
explain below. 

7.2 Temporality and the Crane Model for Quantum 
Gravity 

Crane uses a partition function defined for a triangulated four-manifold. Let 
us denote the partition function by Z(M A ,A, B) =< A\B > M where M 4 
is a four-manifold and A and B are (colored - see the next sentence) three 
dimensional submanifolds in the boundary of M. The partition function is 
constructed by summing over all colorings of the edges of a dual complex 
to this triangulation from a finite set of colors that correspond to certain 
representations of the the quantum group U q {SU{2)) where q is a root of 
unity. The sum is over products of 15J q symbols (natural generalizations of 
the 6 J symbols in angular momentum theory) evaluated with respect to the 
colorings. The specific form of the partition function (here written in the 
case where A and B are empty) is 

Z(M 4 ) = ^- e S A n CT dzm (? (A( ( T))n T ^m~ 1 (A(r))n c 15J g (A(0). 

Here A denotes the labelling function, assigning colors to the faces and 
tetrahedra of M 4 and v — e is the difference of the number of vertices and 
the number of edges in M 4 . Faces are denoted by a, tetrahedra by r and 
4-simplices by (. We refer the reader to jH] for further details. 

In computing Z(M 4 , A, B) =< A\B >m one fixes the choice of coloration 
on the boundary parts A and B. The analog with quantum gravity is that a 
colored three manifold A can be regarded as a three manifold with a choice 
of (combinatorial) metric. The coloring is the combinatorial substitute for 
the metric. In the three manifold case this is quite specifically so, since the 
colors can be regarded as affixed to the edges of the simplices. The color on 
a given edge is interpreted as the generalized distance between the endpoints 
of the edge. Thus < A\B >m is a summation over "all possible metrics" 
on M 4 that can extend the given metrics on A and B. < A\B >m is an 
amplitude for the metric (coloring) on A to evolve in the spacetime M 4 to 
the metric (coloring) on B. 
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The partition function Z(M 4 , A, B) =< A\B > M is a topological in- 
variant of the four manifold M 4 . In particular, if A and B are empty (a 
vacuum- vacuum amplitude), then the Crane- Yetter invariant, Z(M A ), is a 
function of the signature and Euler characteristic of the four-manifold [H]. 
On the mathematical side of the picture this is already significant since it 
provides a new way to express the signature of a four-manifold in terms of 
local combinatorial data. 

From the point of view of a theory of quantum gravity, Z(M 4 , A, B) =< 
A\B >Mi as we have described it so far, is lacking in a notion of time and 
dynamical evolution on the four manifold M 4 . One can think of A and B 
as manifolds at the initial and final times, but we have not yet described a 
notion of time within M 4 itself. 

Crane proposes to introduce time into M 4 and into the partition function 
< A\B >m by labelling certain three dimensional submanifolds of M 4 with 
special grouplike elements from the quantum group U q (SU (2)) and extending 
the partition function to include this labelling. Movement across such a la- 
belled hypersurface is regarded as one tick of the clock. The special grouplike 
elements act on the representations in such a way that the partition func- 
tion can be extended to include the extra labels. Then one has the project 
to understand the new partition function and its relationship with discrete 
dynamics for this model of quantum gravity. 

Lets denote the special grouplike element in the Hopf algebra G = U q (SU (2)) 
by the symbol J. Then, as discussed at the end of the previous section, one 
has that the square of the antipode S : G — ► G is given by the formula 
S 2 (x) = J~ x xJ. This is the tick of the clock. The DOC derivative in the 
quantum group is given by the formula DX = [X, J] = J(S 2 (X) — X). I pro- 
pose to generalize the discrete ordered calculus on the quantum group to a 
discrete ordered calculus on the four manifold M 4 with its hyperthreespaces 
labelled with special grouplikes. This generalised calculus will be a useful 
tool in elucidating the dynamics of Crane's model. Much more work needs 
to be done in this domain. 
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8 Appendix on Iterants 



The primitive idea behind an iterant is a periodic time series or "waveform" 

• • • abababababab 

The elements of the waveform can be any mathematically or empirically 
well-defined objects. We can regard the ordered pairs [a, b] and [b, a] as ab- 
breviations for the waveform or as two points of view about the waveform (a 
first or b first). Call [a, b] an iterant. One has the collection of transformations 
of the form T[a,b] = [ka,k~ l b\ leaving the product ab invariant. This tiny 
model contains the seeds of special relativity, and the iterants contain the 
seeds of general matrix algebra! Since this paper has been a combination of 
discussions of non-commutativity and time series, we include this appendix 
on iterants. A more complete discussion will appear elsewhere. For related 
discussion see p^l UHl H31 [2B 123 EU E21 E] • 

Define products and sums of iterants as follows 

[a, b] [c, d] = [ac, bd] 

and 

[a, b] + [c, d] = [a + c,b + d]. 

The operation of juxtapostion is multiplication while + denotes ordinary 
addition in a category appropriate to these entities. These operations are 
natural with respect to the structural juxtaposition of iterants: 

...abababababab... 

...cdcdcdcdcdcd... 

Structures combine at the points where they correspond. Waveforms combine 
at the times where they correspond. Iterants conmbine in juxtaposition. 

If • denotes any form of binary compositon for the ingredients (a, b,...) of 
iterants, then we can extend • to the iterants themselves by the definition 
[a, b] • [c, d] = [a • c,b • d). In this section we shall first apply this idea to 
Lorentz transformations, and then generalize it to other contexts. 

So, to work: We have 

[t-x,t + x] = [t,t] + [-x,x] = t[l, 1] + x[-l, 1]. 
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Since [1, l][a, b] = [la, lb] = [a, b] and [0, 0][a, 6] = [0, 0], we shall write 

1 = [1,1] 

and 

= [0,0]. 

Let 

a =[-1,1]. 

(7 is a significant iterant that we shall refer to as a polarity. Note that 

GO = 1. 

Note also that 

[£ — x, i + x] = £ + xa. 

Thus the points of spacetime form an algebra analogous to the complex 
numbers whose elements are of the form £ + xa with aa = 1 so that 

(£ + xa){t' + x'a) = tt' + xx' + (tx' + xt')a. 

In the case of the Lorentz transformation it is easy to see the elements of the 
form [k, k" 1 ] translate into elements of the form 

T(V) = [(1 + v)/y/(l-V*), (1 - v)/y/(l-V*)] = [k, AT 1 ]. 

Further analysis shows that v is the relative velocity of the two reference 
frames in the physical context. Multiplication now yields the usual form of 
the Lorentz transform 

T k {t + xa) = T(v){t + xa) 
= (iM 1 -^ 2 ) - va/y/(l-v*))(t + xa) 



(£ - XV) /^J(l-V 2 ) + (X- vt)cr/yJ{l-V 2 ) 

= £' + x'a. 



The algebra that underlies this iterant presentation of special relativity 
is a relative of the complex numbers with a special element a of square one 
rather than minus one (i 2 = —1 in the complex numbers). 
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The appearance of a square root of minus one unfolds naturally from 
iterant considerations. Define the "shift" operator D on iterants by the 
equation 

D[a,b] = [b,a]. 

Sometimes it is convenient to think of D as a delay opeator, since it shifts 
the waveform ...ababab... by one internal time step. Now define 

i[a,b] = aD[a,b] = [-1, l][b,a] = [-b,a]. 

We see at once that 

ii[a, b] = [—a, —b] = [— 1, — l][a, 6] = (— l)[a, &]. 

Thus 

ii = —1. 

This is the traditional construction of the square root of minus one in terms 
of operations on ordered pairs. Here we have described i[a, b] in a new way 
as the superposition of the waveforms a = [—1,1] and D[a,b] where D[a,b] 
is the delay shift of the waveform [a, b]. 

8.1 MATRIX ALGEBRA VIA ITERANTS 

Matrix algebra has some strange wisdom built into its very bones. Consider 
a two dimensional periodic pattern or "waveform." 



. . . abababababababab. . . 
. . . cdcdcdcdcdcdcdcd. . . 
. . . abababababababab. . . 
. . . cdcdcdcdcdcdcdcd. . . 
. . . abababababababab. . . 



(a b \ ( b a \ ( c d \ ( d c \ 
c d )\d c )\a b )\b a) 
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Above are some of the matrices apparent in this array. Compare the matrix 
with the "two dimensional waveform" shown above. A given matrix freezes 
out a way to view the infinite waveform. In order to keep track of this 
patterning, lets write 

[a, d] + [b, c]r) = 



c d 



where 
and 

( 1 \ 



77 



1 



The four matrices that can be framed in the two-dimensional wave form are 
all obtained from the two iterants [a, d] and [b, c] via the delay shift operation 
D[x, y] = [y,x] which we shall denote by an overbar as shown below 



D[x,y] = [x,y] = [y,x\. 

Letting A = [a,d] and B = [b,c], we see that the four matrices seen in the 
grid are 

A + Br], B + Arj, B + Ar], A + Br}. 

The operator 77 has the effect of rotating an iterant by ninety degrees in the 
formal plane. Ordinary matrix multiplication can be written in a concise 
form using the following rules: 

7777 = 1 

vQ = Qv 

where Q is any two element iterant. 

For example, let e = [—1, 1] so that e = — e and ee = [1, 1] = 1. Let 

i = er]. 

Then 

ii = erjerj = elr/r] = e(— e) = — ee = —1. 
We have reconstructed the square root of minus one in the form of the matrix 
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i = eq= [— 1, 1]t7 = 

More generally, we see that 

(A + Br])(C + Dr]) = (AC + BD) + (AD + BC)r] 

writing the 2x2 matrix algebra as a system of hypercomplex numbers. Note 
that 

(A + Brj)(A - Br]) = AA - BB 

The formula on the right corresponds to the determinant of the matrix. Thus 
we define the conjugate of A + Br] by the formula 

A + Br] = A- Br]. 

These patterns generalize to higher dimensional matrix algebra. 

It is worth pointing out the first precursor to the quaternions: This pre- 
cursor is the system 

{±l,±e,±J7,±i}. 

Here ee = 1 = r/r] while % = er] so that ii = —1. The basic operations in 
this algebra are those of epsilon and eta. Eta is the delay shift operator 
that reverses the components of the iterant. Epsilon negates one of the 
components, and leaves the order unchanged. The quaternions arise directly 
from these two operations once we construct an extra square root of minus 
one that commutes with them. Call this extra root of minus one \J — 1. Then 
the quaternions are generated by 

{i = er), j = yf-Le, k = yZ-Lr)} 

with 

i 2 = j 2 = k 2 = ijk = — 1. 

The "right" way to generate the quaternions is to start at the bottom iterant 
level with boolean values of and 1 and the operation EXOR (exclusive or). 
Build iterants on this, and matrix algebra from these iterants. This gives the 
square root of negation. Now take pairs of values from this new algebra and 
build 2x2 matrices again. The coefficients include square roots of negation 
that commute with constructions at the next level and so quaternions appear 
in the third level of this hierarchy. 



( : »' 
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8.2 Matrix Algebra in General 

Construction of matrix algebra in general proceeds as follows. Let M be an 
nxn matrix over a ring R. Let M = (m^) denote the matrix entries. Let n be 
an element of the symmetric group S n so that 7Ti, 7T2, • • • , 7r n is a permuation 
of 1, 2, • • • , n. Let v = (i>i, v 2, • • • , v n ) denote a vector with these components. 
Let A(v) denote the diagonal matrix whose i — th diagonal entry is i>j. Let 
v n = (v wi , ■ ■ ■ ,v 7Tn ). Let A 7r (t>) = A(v w ). Let A denote any diagonal matrix 
and A n denote the corresponding permuted diagonal matrix as just described. 
Let [n] denote the permutation matrix obtained by taking the i — th row of 
[n] to be the 7Tj — th row of the identity matrix. Note that [7r]A = A 71 " [71"]. For 
each element n of S n define the vector v (M, 71) = (m l7ri , • • • , m n7Tn ) and the 
diagonal matrix A[M] n = A(v(M,ir)). 

Theorem. M = (l/(n - l)!)E 7r65 „A[M] 7r [7r]. 

The proof of this theorem is omitted here. Note that the theorem expresses 
any square matrix as a sum of products of diagonal matrices and permu- 
tation matrices. Diagonal matrices add and multiply by adding and multi- 
plying their corresponding entries. They are acted upon by permutations as 
described above. This means that any matrix algebra can be embedded in 
an algebra that has the structure of a group ring of the permutation group 
with coefficients A in an algebra (here the diagonal matrices) that are acted 
upon by the permutation group, and following the rule [vr]A = A 71 " [71-]. This 
is a full generalization of the case n = 2 described in the last section. 

It is amusing to note that this theorem tells us that up to the factor of 
l/(n — 1)! a unitary matrix that has unit complex numbers as its entries is 
a sum of simpler unitary transformations factored into diagonal and permu- 
tation matrices. In quantum computing parlance, such a unitary matrix is 
a sum of products of phase gates and products of swap gates (forming the 
permutations). 

A reason for discussing these formulations of matrix algebra in the present 
context is that one sees that matrix algebra is generated by the simple op- 
erations of juxtaposed addition and multiplication, and by the use of per- 
mutations as operators. These are unavoidable discrete elements, and so the 
operations of matrix algebra can be motivated on the basis of discrete phys- 
ical ideas and non-commutativity. The richness of continuum formulations, 
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infinite matrix algebra, and symmetry grows naturally out of finite matrix 
algebra and hence out of the discrete. 



9 Philosophical Appendix 

The purpose of this appendix is to point to a way of thinking about the 
relationship of mathematics, physics, persons, and observations that under- 
lies the approach taken in this paper. We began constructions motivating 
non-commutativity by considering sequences of actions ■ ■ ■ DCBA written 
from right to left so that they could be applied to an actant X in the order 
• • • DCBAX = ■■■ (D(C(B(AX))) • • • . The sequence of events A, B, C, D, ■ ■ ■ 
was conceptualized as a temporal order, with the events themselves happen- 
ing at levels or frames of successive "space". There is no ambient coordi- 
nate space, nor is there any continuum of time. All that is given is the 
possibility of structure at any given moment, and the possibility of distin- 
guishing structures from one moment to the next. In this light the formula 
DX = [X, J] = XJ — JX = J(X' — X) connotes a symbolic representation 
of the measurment of a difference across one time interval, nothing more. 
In other words DX represents a difference taken across a background differ- 
ence (the time step). Once the pandora's box of measuring such differences 
has been opened, we are subject to the multiplicities of forms of difference 
Vif X = [X, K] , their non-commutativity among themselves, the notion of a 
flat background that has the formal appearance of quantum mechanics, the 
emergence of abstract curvature and formal gauge fields. All this occurs in 
these calculi of differences prior to the emergence of differential geometry or 
topology or even the notion of linear superposition of states (so important 
to quantum mechanics). Note that in this algebraic patterning each algebra 
element X is an actant (can be acted upon) and an actor (via the operator 
Vx). In Lie algebras, this is the relationship between the algebra and its 
adjoint representation that makes each element of the algebra into a repre- 
sentor for that algebra by exactly the formula aajA(X) = [A,X] = —Va(X) 
that we have identified as a formal difference or derivative, a generator for a 
calculus of differences. 

The precursor and conceptual background of our particular formalism 
is therefore the concept of discrimination, the idea of a distinction. A key 
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work in relation to that concept is the book "Laws of Form" by G. Spencer- 
Brown in which is set out a calculus of distinction of maximal simplicity 
and generality. In that calculus a mark (denoted here by a bracket < >) 
represents a distinction and is seen to be a distinction between inside and 
outside. In this elemental mathematics there is no distinction except the 
one that we draw between the mathematician and the operator in the for- 
mal system as sign/symbol/interpretant. This gives full responsibility to the 
mathematician to draw the boundaries between the formal system as phys- 
ical interaction and the formal system as symbolic entity and the formal 
system as Platonic conceptual form. In making a mathematics of distinc- 
tion, the mathematician tells a story to himself /herself about the creation 
of a world. Spencer-Brown's iconic mathematics can be extended to contact 
any mathematics, and when this happens that mathematics is transformed 
into a personal creation of the mathematician who uses it. In a similar (but 
to a mathematician) darker way, the physicist is intimately bound to the 
physical reality that he studies. 

We could have begun this paper with the the Spencer-Brown mark as 
bracket: < > . This empty bracket is seen to make a distinction between 
inside and outside. In order for that to occur the bracket has to become 
a process in the perception of someone. It has to leave whatever objective 
existence or potentiality it has alone (all one) and become the locus or nexus 
of an idea in a perceiving mind. As such it is stabilized by that percep- 
tion/creation and becomes really a solution to{< >}=< > where the 
curly bracket (the form of perception) is in the first place identical to the 
mark < >, and then distinguished from it by the act of distinguishing world 
and perceiver. It is within this cleft of the infinite recursive and the finite 

< >={< >} = {{< >}} = {{{< >}}} = ■•• = {{{{{{••■}}}}}} 

that the objectivity of mathematics/physics (they are not different in the 
cleft) arises. All the rest of mathematics or calculus of brackets needs come 
forth for the observer in the same way. Through that interaction there is 
the possibility of a deep dialogue of many levels, a dialogue where it is seen 
that mathematics and physics develop in parallel, each describing the same 
boundary from opposite sides. That boundary is the imaginary boundary 
between the inner and outer worlds of an individual. 
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